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Abstract 

Antiferromagnetic and charge ordered Hartree-Fock solutions of the one-band Hubbard model 
with on-site and nearest-neighbor Coulomb repulsions are exactly mapped onto an auxiliary local 
Kohn-Sham (KS) problem within a density- functional theory. The mapping provides a new insight 
into the interpretation of the KS equations, (i) With an appropriate choice of the basic variable, 
there is a universal form of the KS potential, which is applicable both for the antiferromagnetic and 
the charge ordered solutions, (ii) The Kohn-Sham and Hartree-Fock eigenvalues are interconnected 
by a scaling transformation, (iii) The band-gap problem is attributed to the derivative discontinuity 
of the basic variable as the function of the electron number, rather than a finite discontinuity of the 
KS potential, (iv) It is argued that the conductivity gap and the energies of spin-wave excitations 
can be entirely defined by the parameters of the ground state and the KS eigenvalues. 

PACS numbers: 71.15.Mb, 71.45. Gm, 71.45.Lr, 75.30.Fv 
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I. INTRODUCTION 



Practical success of the Kohn-Sham (KS) scheme of density-functional theory (DFT)^ 
in theoretical description of the ground-state properties of molecules and solids lies in the 
idea of substitution of the real interacting many-electron system with an auxiliary system 
of noninteracting particles moving in an effective local potential. Similar ideas can be found 
recently in the dynamic mean-field theory of strongly correlated systems,'^ which is based 
on a mapping of lattice models (e.g., the Hubbard model or the Kondo model) onto a lo- 
cal single-impurity problem. This mapping is well elaborated only for a limited number of 
examples like the homogeneous electron gas in the density-functional theory [the so-called 
local-spin-density approximation (LSDA)], or in the limit of infinite lattice coordination in 
the dynamic mean-field theory. However, many general aspects of such a local noninteracting 
representation underlying the KS scheme are not well understood. One of the most contro- 
versial issues is the physical meaning (if any) of the KS eigenvalues, and their relationship 
with true physical one-particle excitations. 

The major problem, which was intensively discussed in the 1980s ,^i^i^i'''i^i^i^°i^'^i'^^i'^^ is the 
description of fundamental conductivity gaps in semiconductors and insulators. It is gen- 
erally accepted that the KS one-electron orbital energies have no precise physical meaning, 
except for the highest occupied eigenvalue, which is equal to the chemical potential (minus 
the work function). Formally, the conductivity gap is not the ground-state property and 
includes the excited-state effects. The LSDA systematically underestimate the band gaps. 
The band-gap problem is typically ascribed to the finite discontinuity of the KS potential 
as the function of the electron number, meaning that for an extra electron, the KS potential 
may have an additional constant compared to the KS potential designed for the ground 
state. In other words, this is equivalent to the statement that the band-gap problem for 
semiconductors and insulators may be resolved by a rigid upward shift of the conduction 
bands in the KS picture for the ground state. This procedure is frequently referred to as a 
" scissor operator" . 

None of the standard density functionals, such as the LSDA or generalized gradient 
expansions show such a discontinuity. The origin and the magnitude of this discontinuity 
is still an open question. Several authors argued that the discontinuity can be large, and 
the single KS potential (even an exact one) designed for the ground state does not give 
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the correct gap .^i^i'^'i^ In accordance to their opinion, the solution for the band-gap problem 
should be sought in the nonlocality and the energy- dependence of the electron self-energy, 
which is beyond the standard KS formulation, and can be obtained, for example, in the 
framework of the GW approximation /i^i^i^^ On the contrary, Gunnarsson and Schonhammer 
suggested that the actual discontinuity of the exact KS potential may be small, and the main 
problem lies in the LSDA.— Their original arguments were based on a numerical analysis of 
a Hubbard-type model with the on-site Coulomb repulsion only. However, later they also 
argued that the intersite Coulomb interaction can produce a large discontinuity of the KS 
potential.— The disputed issues have been addressed again recently in Refs. QandH usmg 
two different viewpoints on the formulation of the density functional-theory on a discrete 
lattice. 

In this work we discuss several general aspects of the density-functional theory for the one- 
band extended Hubbard model. We focus our attention on the analysis of antiferromagnetic 
(AF) and charge ordered insulating states, the description of which was a traditionally diffi- 
cult subject for ab initio DFT methods . ^^'^^'^^ From the very beginning we restrict ourselves 
by the framework of the (static) mean-filed (MF) approximation. Physically, this means 
that we take into the consideration only the (classical) Coulomb and exchange interactions 
among the electrons, and neglect the correlation effects. In this sense our density functional 
theory for the Hubbard model is certainly less advanced rather than similar theories exist- 
ing in the literature and dealing with a formally exact problem numerically.— On the 
other hand, we hope that such a simplified mean-filed analysis may be rather instructive, 
because we will be able to treat the problem analytically and give a clear picture for several 
disputed moments. Our primary goal will be to show how the same mean-field problem can 
be formulated in two different forms, one of which is the conventional Hartree-Fock (HF) 
theory, where the one-electron orbitals have a clear physical meaning based on Koopmans' 
theorem, and the other one is the KS scheme in the DFT. The main point, which we would 
like to emphasize in this work is that the meaning of the KS band structure may not be so 
formal, as it is currently understood. The relationship between the KS eigenvalues and the 
true one-particle excitations may be more fundamental, and is not necessarily based only 
on the discontinuity property of the exact KS potential. We do not discuss here the relative 



stability of different solutions o: 
found elsewhere (see, e.g., Refs 



: the extended Hubbard model. Such an analysis can be 
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The paper is organized as follows: the total energy functional is defined in SecjlTl The KS 
potential is derived in Sec JIlII by using the optimized effective potential (OEP) method.— 
The connection between the OEP method, the density functional theory and the unrestricted 
Hartree-Fock approach is discussed in Sections II V I and IVI In Sec iVll we discuss the connection 
between the conductivity gap, the energies of spin-wave excitations and the KS one-electron 
energies. Brief summary will be given in Sec lVIII 

II. TOTAL ENERGY 

The one-band extended Hubbard Hamiltonian with on-site {U) and nearest-neighbor (V) 
Coulomb repulsions is defined as 

^ = E (^^^ + ^^^) + ^ E ^I^i + ^ E (1) 

where n1=cl'^'5l and c^^ (c"") creates (annihilates) an electron on the site i with spin a=] or 
[. A nearest-neighbor pair is denoted by {ij). The kinetic hopping between nearest-neighbor 
orbitals is denoted by t. We use the convention, in accordance with which capital letters are 
reserved for absolute values of parameters of Coulomb interactions and the KS potential, 
whereas small letters stand for the same parameters in units of 2t (for example, u=U/2t, 
etc.). 

Let us consider the mean-field solution of the extended Hubbard model, where the ground- 
state wavefunction of the many-electron system is approximated by a single Slater deter- 
minant of one-electron orthonormal orbitals {ipm}- Then, the total energy is given 
by 

SuF[rj = {s[rj\n\s[rj), (2) 

and can be expressed in terms of the one-electron density matrix 

p''[^,j] = J2^rj-^^'^\rj (3) 

m 

as £^MF=^+^H+^^x (see, e.g., Ref. Q), where T is the (noninteracting) kinetic energy: 

r = -tY,ipy,J]+P''[J,^]), (4) 
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is the Hartree (Coulomb) energy: 
and £x is the exchange (Fock) energy: 

^^ = -^11 ^p""^'^ ^])' - ^ E P'^[^'^>'^b-, (6) 

j(T (ij) 

u 

Both and include the self-interaction. 

Our goal is to simulate a spin or /and charge ordering on the bipartite one-dimensional 
(ID), two-dimensional (2D), simple cubic (SC) or base-centered cubic (BCC) lattice, for 
which the spin or/and charge densities are periodically distributed between two types of 
atoms v= 1, 2 resulting in the doubling of the unit cell, as it is shown in Fig^ Thus, 
the site-indices are specified by the position of the doubled unit cell R and the type of the 
site: ^=(Rz/). p" is the periodic object: p'^[RV, Rz/]=p'^[(R' — R)z/',Oz/], implying that the 
one-electron orbitals in Eq.© are the Bloch states: i.e., m=(k£) and c^^|V'k£)=e*''^c^y|'?/'^^), 
where k is a vector in the first Brillouin zone of the doubled unit cell, and I is the remaining 
(band) index. In accordance with the definition (^, the range of interactions in the Hamil- 
tonian Ti is restricted by the nearest neighbors. If we combine this fact with the symmetry 
properties of the density matrix, i.e., the periodicity and the invariance with respect to the 
local symmetry operations for the lattices shown in Fig^ it is rather easy to see that the 
total energy Q will be entirely represented by the elements of the density matrix with 
R=R'= 0, which can be formally composed into the 2x2 matrix: 

^^1 p'^[01,01] p'^[01,02] 
\^p'^[02,01] p'"[02,02] 

Since the inversion is one of the local symmetry operations, p°" may be taken in the real 
symmetric form. 

In the following we will focus out attention on the three possible ordered states: 
the antiferromagnetic state near to half filling, 
the charge-ordered nonmagnetic (CON) state near to half filling, 
the charge-ordered ferromagnetic (COF) state near to quarter filling. 

The number of electrons in the doubled unit cell is specified as ne=no±?7, where hq is the 
integer electron number at half/quarter filling (no=2 for the AF and CON state, and no=l for 
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the COF state) and t] describes a small deviation away from half /quarter filling. Then, the 
elements of the density matrix for these three cases obey the series of symmetry constraints 
listed in Table El reducing the number of independent elements in p^'^, for a given rie, to the 
expectation values of only two pseudospins calculated for the majority-spin channel: 

Tx = Tr(pfax), = Tr(p^aJ, (8) 

where ax and are 2x2 Pauli matrices. In the AF state, is equal to the local magnetic 
moment. If we define the number of the doubled unit cells as A^, and the total energy 
per doubled unit cell as Emf=S/N=T+Eii+Ex, we obtain from Eqs.ijH)-® ^^e following 
expressions for the kinetic energy (T), the Hartree energy (-Eh), and the exchange energy 
(Ex) per doubled unit cell, correspondingly 

T[t^] = -noztr^, (9) 

Mr.] = \nl{U + zV) + \{U - zV){t. + ri)^ (10) 



and 

1 { Tl 

Ex[ry,,T^] = --uqU 



\n,zVrl (11) 



,^0 

where z is the coordination number. In Eq. ()10|) . is — r^, and for the AF, CON and 
COF states, respectively. The prefactor hq in Eqs.® and (fTT|l counts the spin contributions. 
Due to the symmetry properties of the density matrix (Table E}, we have to equivalent 
contributions to T and Ex coming from the majority-spin and minority-spin channels in the 
case of AF and CON ordering. This is reflected in the prefactor no=2, which coincides with 
the number of electrons at half filling. In the case of COF ordering, the minoriti-spin states 
are unoccupied. Therefore, we have to count only the majority-spin channel and take no=l, 
which again concides with the number of electrons at quarter filling. 



III. OPTIMIZED EFFECTIVE POTENTIAL 



There are two possible ways of searching the ground state for the mean-field total energy 
functional (j2I)- One is the standard Hartree- Fock approach, which will be discussed in 
SecEl Here we start with an alternative provided by the KS scheme and assume that the 
one-electron orbitals {V'k^} can be obtained as UeXN lowest energy eigenf unctions of the 
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following Hamiltonians: 

(RIR'2) R 

As in Eq.(0), the notation (RIR'2) means that the summation is restricted by the nearest 
neighbors, which are now of different types (see Fig^. If r is the relative position of two 
different atoms within one unit cell, the notation (RIR'2) is equivalent to the condition 
|R' + r — R| = |r|. Thus, ?i]^g retains the kinetic part of by the original Hamiltonian Ti.. 
The interacting part is substituted by the effective one-electron KS potentials {B'^}, which 
are assumed to be site-diagonal ("local" with respect to the site- indices), do not depend 
on R (as it is requested by the periodicity of the antiferromagnetic and charge-ordered 
states), but may depend on the type of the atom u. {B^} can be found in the spirit of 
optimized effective potential metho d^^i^^ as the solution of the variational problem (for the 
given electron number ne): 

d 



dBz 



TI 



0, (13) 



provided that {B^} obey the series of constraints jF^fji?^}], which guarantee the required 
symmetry of the density matrix (Table |l}. Thus, {-B|^} can be interpreted as the exact KS 
potentials for exchange alone.— The Hamiltonian (fT^ can be diagonalized in the basis of 
Bloch orbit als 

1 -R\y (14) 

c^t(k); ViVV {e^'-'c^, ' 




(15) 



The eigenfunctions of Ti^g are 

/ ~ \-sin^9^ cost?^ ) \^,\m 
where |) denotes the vacuum state with no electrons in it. 

We note that due to the symmetry constraints (Table H]), we have to consider only the 
case —B\=B\=B. Then, the eigenvalues of ifj^g are 

4(k) = ±v/52 + 4t2/'(k) (16) 

and the angle §\. in Eq.()15|) is given by 

tg^T ^ fi^) (17) 
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where 

/(k) = i Yl cosk(r-R) 

R:|r-R| = |r| 

is the structure factor. [The exphcit from of /(k) for different lattices is shown in FigC]. 
Then, putting the wavefunctions (fT^ in equations ((Tj) and (jH)), we obtain after some 
algebra: 

r^ = ^i: 4==T0N - elm (18) 



zN 



-E 



-i)b 



-.Q[sF-slik)], 



where £=±1, and the Fermi energy e^? is determined by the condition: 



(19) 



(20) 



Then, it is rather straightforward to show for any fixed electron number rie that 



db 



2b dr, 
z db 



(21) 



Therefore, Tx and are not independent variables and the variational procedure of the OEP 
approach is equivalent to a DFT in which participates as the basic variable. Indeed, using 



Eqs.® and (EH) we obtain {dT[T^]/di 



-uqB. As it was discussed earlier, the prefactor 



no counts the spin channels involved into the problem. Then, variational principle leads 
to the standard definition for the KS potential conjugated with the basic variable t^: 



no 



9r, 



(22) 



where the minus sign is related with our definition of and B: for cr=t, the solution with 
rz=p^ [01, 01] [02, 02] >0 corresponds to the potential bI=-BI=-B<0. The exchange 
part in Eq. (P^ is given by 



Bx^ 

no 



1 dEx[T^,T^ 



2 



Vr^b 



(23) 



[see Eqs.dll]) and The Coulomb part B^={-l/no){dEYi[T^]/dT^)\n, is 0, (U-zV)t, 

and \{U-zV)t^ for the AF, CON and COF states, respectively [see Eq.jTUl)]. Then, Eq. 
can be transformed to the self-consistent equation for the KS potential 6=6h+&x: 

1 



b = -mhf1"z 



VT^b, 



(24) 



where mhf is u, 2zv — u and zv for the AF, CON and COF states, respectively. Eq. 
further yields: 

& = ^mks[i"x]^z, (25) 
where the effective coupling uks[tx] is given by 

"KslrJ = (26) 

Putting Eq. ()25|l in Eq.((2SI), one can find for the self-consistent exchange potential hx'- 

1 



(27) 



1 + VT^ 

The local KS potential has the same form [Eq. (j25p ] for different ordered states, although 
the driving force for the ordering may be different: u in the AF state and v in the CO 
states. The local form of the KS potential for the CO states has one important implication 
which we would like to emphasize here. It is well established that Hubbard-type models 
with solely on-site Coulomb interactions treated on the level of mean-field theories may 

n 

exhibit the charge-ordered solutions (see, e.g., Refl23). The present work suggests that such 
"on-site Coulomb interactions" may be interpreted in the Kohn-Sham sense, as parameters 
of an auxiliary local KS problem, whose real physical meaning may be different and implies 
the contributions of longer-range Coulomb interactions. 

uks has the common renormalization factor (l+fTx)^^ which is driven by the intersite 
Coulomb interaction v and depends on through t^. The dependence Ty^r^ is shown 
in Fig 121 for several lattices and ne=no. We note: (i) for a given r^, is larger when z 
is smaller. Thus, the renormalization is the most important for ID systems, (ii) In the 
framework of our model, Tz=0 may be regarded as a homogeneous solution. Since Tx is 
the decreasing function of r^, we have uks[0]<mks[''"z] meaning that the parameter uks being 
once evaluated in the homogeneous limit and employed afterwards for the whole range 
of Tz is generally " overscreened" . It is therefore clear that an attempt to substitute the 
real effective coupling uks[tz] with its homogeneous counterpart uks[0] in the KS potential 
will generally underestimate the tendencies towards the antiferromagnetism and the charge 
ordering, as it is well established for the LSDAj^ The analytical estimate for Tz=0 and 
ne=no is rx[0] = (2/7r)'°S2 2 (except the SC lattice). At the opposite end rz=l, corresponding 
to the limit nKs[l]^C)0, we have rx^l/26-^l/nKs[l]5 see Eqs. (fTH|) and (j^ . Then, from 
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Eq.(j2ni) we obtain MKs[l]='UHF(l+'y/wKs[l]) "^7 whence it follows that 

mks[1] = mhf - V. 



(28) 



IV. THE LEVY CONSTRAINED-SEARCH FORMULATION 

A very elegant formulation of the density-functional theory was introduced by Levy.— For 
our purposes it is read as follows: if is the basic variable, the total energy minimization 
(for a given electron number) can be performed in two steps: 

8 = min \ min (^'|H|\E') 

Tz [|*)^Tz 

The inner minimization is constrained to all many-electron wavefunctions that give r^, 
while the outer minimization releases the constraint by searching all r^. In the mean-filed 
approximation, the search of all possible many-electron wavefunctions is restricted within 
the subspace of Slater determinants Thus, is additionally constrained, and 

we have S<Smf- The purpose of this section is to show that in all other respects, the 
OEP approach employed in the previous section is equivalent to the Levy constrained- 
search formulation (j^^ . The proof is simple and can be carried out along the general 
line taken in the standard DFT formulation.^ Taking into account the explicit form of 
('^'[V'm] I'^l'S'i^m])' Eqs.(P)- (fTT|) . the inner step in Eq. ipi?]) is equivalent to the search of the 
one-electron orbitals {ip!^} which maximize the expectation value of the pseudospin for 
the given r^, or equivalently, minimize the (noninteracting) kinetic energy T for the given 
Tz. Using the definitions of T and in terms of the one-particle orbitals [see Eqs. ©,(111), 
((Tj) and ©], this gives the functional 

m (RIR'2) m R mm' 

(30) 

and the minimization condition 6Q['iplj]/6{ip}^\=0. (The Lagrange multipliers e]^^, stand for 
the orthonormalization condition between IV^^'))- These result in the equation 

h'\i^L)=Y.^m'\^l')^ (31) 

m' 

where 

— ~^ ('^Rl'^R'2 + '^R'2'^Rl) ^ ('^Rl'^Rl ~ '^R2'^R2) • (^2) 

(RIR'2) R 
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Similar equations can be obtained for cr=| by using the symmetry properties of the density 
matrix from Tabled The Hamiltonian (j^^ has the same form as the KS Hamiltonian (fT^ . 
(Note also that since is Hermitian, sl^^, can be taken in the diagonal form). If is 
the expectation value of the pseudospin corresponding to the ground state, the Lagrange 
multiplier is equal to the self-consistent potential B. Thus, the inner minimization 
of the Levy constrained-search formulation is equivalent to the diagonalization of the KS 
Hamiltonian in the OEP approach. 

The conclusion is very natural and, in principle, directly follows from the definition of the 
KS scheme.-^ Indeed, the mean-field energy £^mf is determined by only two parameters, 
and Tz. Both of them are requested to come out exactly of the one-electron KS orbitals {"^A^}- 
Tz as the basic variable, and Tx as the only variable which determines the noninteracting 
kinetic energy. 

V. UNRESTRICTED HARTREE-FOCK APPROACH 

In order to clarify the origin of the renormalization factor (l+fTx)^^ in the KS potential 
fl25|) . let us now turn to the Hartree-Fock solution of the extended Hubbard model, where the 
basic variables are the one-electron orbitals {ipm}'- '^HFlV'rn)='^^MF/'5(V^ml- Taking into ac- 
count the periodicity of the ordered states, this gives the Hamiltonian (the uniform potential 
shift is dropped): 

^HF = -^(1 + ^^x) Yl (^Ri^R'2 + ^rAi) - ^Uuft, J2 (^RiCri " CR2CR2) (33) 

(RIR'2) R 

with the parameters of site-diagonal interactions ?7hf defined earlier. The prefactor {1+vt^) 
in the off-diagonal part of the HF Hamiltonian comes from the nonlocality of the Fock 
exchange potential, which is combined with the kinetic term. Since the KS potential is re- 
quested to be local, the same effect enters the KS formulation through the renormalization 
of the site-diagonal interactions, as it is indeed manifested in Eq. fl26|) . Therefore, the termi- 
nology "the local potential" and "the local approximation" used in the density-functional 
theory does not necessarily mean that the nonlocal exchange interactions are totally ne- 
glected. They may contribute to the renormalization of the local interactions. 

In the framework of the one-band extended Hubbard model, we obtain the trivial connec- 
tion between the KS and HF Hamiltonians: 'Hl^g={l+VTx)~^Ti.lj^p, meaning that they have 
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the same set of eigenfunctions, but different eigenvalues. The KS band structure is simply 
the scaled HF one. Since Tx>0 [see Eq. fllSj) ]. the scaling is in fact a contraction, unless v=0. 
Although such a trivial scaling scenario is a crude simplification of the situation realized in 
realistic calculations, it does explain the main tendency: not only the band gaps, but also 
the bandwidths appear to be smaller in the KS picture.— This is definitely different from 
the rigid upward shift of the conduction bands, underlying the "scissor operator", which 
affects only the gap width. 

VI. ONE-PARTICLE EXCITATIONS 

In this section we establish several useful connections between Kohn-Sham eigenvalues 
and true one-particle excitation energies for the extended Hubbard Hamiltonian in the mean- 
filed approach. First, we will considered the conductivity gap as an example of excitations 
for which the total number of particles is not conserved. We will show that although in 
this case the physical gap can be obtained from the total energy functional and is entirely 
defined by the parameters of the ground state, the one-electron KS band structure plays 
a formally fictitious role, in the sense that the KS gap is generally different from the true 
conductivity gap. Second, we will turn to the analysis of virtual spin-wave excitations in 
the antiferromagnetic insulating state at half filling and argue that these excitations can be 
interpreted directly in terms of the KS eigenvalues. 

A. Conductivity gap 

If ne=no, the KS scheme yields the insulating solution. The band gap AKs=f^KsTz is direct 
and opens in the points {ko} for which /(ko)=0 and '*9k(,=0 [see Eqs. (fTBj) and (fTTj) ]. In 
order to calculate the physical conductivity gap, we first take the standard procedure^^ and 
employ the fact that two Hamiltonians, Tij^g and T^Hp, have the same set of eigenfunctions. 
Then, the conductivity gap is given by 

A = Aks + (< JT^^F - nls\^l+) - (< J^HF - ^Ksl<-)- (34) 

Taking this as an exercise and using Eq. ()15|) for '&l^^=0, it is very easy to verify that 
A=f/HFTz=AHF; which is nothing but the consequence of well-known Koopmans' theorem 
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for the eigenvalues of the HF Hamiltonian ()33|) . (Note, that the correction to Koopmans' 
theorem due to the relaxation of HF levels vanishes for extensive systems when N-^oo). 

The central question is how the physical gap is related with the ground-state properties, 
and whether it can be entirely described by the parameters of the ground state and the 
KS eigenvalues. As it was discussed in the introduction, the most prominent scenario is 
based on the idea of discontinuity of the KS potential. In the case of finite species the 
similar problem was investigated for the chemical hardness (half of the band gap in the limit 
A^-^oo). It was argued that the exact KS potential exhibits a constant jump (discontinuity) 
when the number of electrons changes through an integer value, and the chemical hardness 
can be found as the difference of highest occupied KS eigenvalues calculated for the initial 



and the final electron configurations (see, e.g., Refs. 33 |27|). The constant jump of the 
potential, however, is not related with the ground-state properties. Similar arguments can 
be applied for infinite systems, where the band gap can be expressed through the change of 
the chemical potential.— In the context of the present work, this would be equivalent to the 
claim that the uniform part of the KS potential (C), which was irrelevant to the search of 
the ground state and dropped out the analysis in Sections IIIHlV| should exhibit the finite 
jump across the band gap: C|„(,+^— C|„q_^=Ahf— Aks, and the band-gap problem might be 
resolved via the scissor-type correction. 

Below we present a different view on the band-gap problem. To begin with we note that 
the quantities t^, and b, which do determine the ground state, exhibit no finite disconti- 
nuity as the function of the electron number (FiglH)). The analysis presented in preceding 
sections clearly shows that the physical gap Ahf can be expressed in terms of these charac- 
teristics of the ground state by using another, scaling transformation AHF=(l+'yTx) Aks-~ 
Thus, there is a unique connection between Ahf and the ground-state properties. In order 
to gain more insight into the problem, we argue below that the physical gap Ahf can be 
naturally derived from the mean-field total energy functional Smf['Tx,'Tz], not exploiting the 
idea of the finite discontinuity of the exact KS potential. Instead, we argue that the band- 
gap problem is related with the derivative discontinuity of the basic variable across the 
band gap. We start with the general definition: 

A = ^mf(^oA^ + 1) + ^MF(noA^ - 1) - 2SMF{noN). (35) 
As in Eq. fj34|) . we assume that the addition or the removal of an electron does not destroy 
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the symmetry of the ground state. (For example, if the ground state of the uqN electron 
system is antiferromagnetic, we assume that the system of uqN^I electrons will remain in 
the antiferromagnetic state, though this requirement may formally contradict to the Nagaoka 
theorem^). Then, in the limit N^oo, Eg . (j35j) becomes 




drip 



dEupiT^, Tz] 



dup 



(36) 



where ri=l/N. Taking into account the explicit form of and r^, Eqs. (jl8p and (jl9|) . we 
will have two contributions to {dEMpiT^, Tz]/ dne)\no±ri- due to the change of the the 

band filling with the fixed potential b (the derivative of the 0-function), and the other one 
is due to the ne-dependence of the self-consistent potential itself. Therefore, we can write: 



drip 



dEMF[r^,rz] 



no±ri 



dup 



+ 



dEMF[r^,r^] 



db 



db 



no 



dup 



no±ri 



However, in the ground state it holds {dEMF['Tx,'Tz]/db)\no=0 [see Eq. (fT^ ]. and we find: 



A 



lim 



5^MFkx,7-z] 



dup 



dEMF[r^,r2] 



b,no+r] 



dup 



(37) 



Then, we have the following identity: (<96[£i;'-£^(k)]/9ne)|6.„(,±,,=5[£F|fe,„o±^-£:±(k)]((9£:i;'/5ne)|b,no±r,- 
In the limit rj—^O, the 5-function will pick up only those points in the Brillouin zone which 
correspond to the band-gap edges, i.e., k G {ko}. Therefore, we obtain from Eq. pn|) : 



lim > 



9e[£p-4(k)] 



dup 



lim > 



dup 



N 
no' 



b,no±ri 



b,no±r] kg{ko} 

This immediately yields: lim^_>o(f^''"x/c^^e)|6,no±r?=0 [see Eq.lfTHj) and note that /(ko)=0], 
meaning that there is no derivative discontinuity of the kinetic energy (jH)) across the band 
gap caused by the change of the band filling with the fixed KS potential. (Do not mix it 
with the behavior of Tx in Fig|21 which shows the total effect, including the change of the 
KS potential b). 

On the contrary, the derivative of the basic variable shows the discontinuity: 
rio lim^^o(<9Tz/9ne)|6,„o±r?==Fl- Then, putting these results in Eq. ljHTj) and using the ex- 
plicit form of Euf['^x,'^z\, Eqs.(P|)-(|lip. we again arrive at Koopmans' theorem for the HF 
gap: A=?7HFrz=AHF- 

Thus, although the KS gap is generally different from the true conductivity gap, Ahf in 
our model, the latter can still be derived from the total energy functional i?MF[Tx,Tz] using 
the definition ^7\. 
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B. Virtual spin-wave excitations 



Let us consider the energies of spin-wave excitations in the antiferromagnetic insulating 
state at half filling. 

For these purposes it is convenient to work with the local representation for the one- 
electron potential provided by the KS scheme. Then, the total energy changes caused 
by small nonuniform rotations of the spin magnetic moments near the (antiferromagnetic) 
equilibrium, i.e., the processes relevant to the spin- wave excitations, can be exactly mapped 
onto the Heisenberg-type model^Si^i*^ 



where ej is the direction of the spin magnetic moment at the site i. The parameters of the 
model can be expressed in terms of the one-electron Green function Glj^{e) constructed for 
the KS Hamiltonian ^ as^ii^ 



where Al^=Bl—Bl is the local exchange splitting at the site ^=(Rz/). The mapping is 
based on so-called local force theorem. One can argue rather generally that if the spin- 
magnetization density (ir^ in our model for the antiferromagnetic state) is the basic variable 
in a (spin)-dinsity-functional theory, and is requested to come out exactly of the one-electron 
KS orbitals, such a mapping onto the Heisenberg-type model should be exact.— 

For two nonequivalent sites of our model we have: A^=—A^=UksTz- Thus, the 
strength of the interatomic exchange couplings Jy is directly related with the magnitude of 
the KS gap AKs=t^KsT-z- 

Let us consider the large-uxs limit. In the lowest order of the 1/uks expansion, the 
interatomic magnetic interactions are restricted by the nearest neighbors (ij). Corresponding 
elements of the Green function are given by G|.^.^(£:)=G'|^.^(£:)^— t/(£:^ — |f/^gr^^). Putting 
them in Eq. ljHIH) and noting that in the large-MKs limit rz=l and Uks=U—V [see Eq.(j2Hl)], we 
arrive at the well-known result for the effective magnetic coupling between nearest neighbors 
in the extended Hubbard model>i2*2a»2ii^ j^..^=t2 /(^U-V), or J(^y)=t^ /Aks- [Note, that the 
sign of J(jj) depends on the convention used for the directions of the spin magnetic moments 
in Eq.()38|). If we refer to the antiferromagnetic ground state, ■ ej=— 1, and therefore 




(38) 




(39) 



15 



J(jj)>0]. Thus, the interatomic magnetic interactions can be regarded as virtual hoppings 
between effective KS levels, and in this sense the KS band structure has the clear physical 
meaning. 



Finally, we would like to mention one rather interesting implication of two examples 
discussed in this section to the theory of strongly-correlated materials. The fact that the 
intersite Coulomb interaction V " renormalizes" primarily the exchange coupling J^jj) in the 
antiferromagnetic state, but not the conductivity gap was well-known from earlier studies 
of the extended Hubbard model. Such an apparently different influence of the intersite 
Coulomb interaction V on different physical quantities has led to some ambiguity in the de- 
termination of the physical parameter, which controls the strength of correlation effects in 
solids (see, e.g., discussions in Refs I21I ). The present work suggests that these two proper- 
ties are not necessarily independent. In the mean-field approach, such a behavior means that 
that there are (at least) two different representations for the one-particle spectrum, which 
correspond to the same ground state. One representation is essentially nonlocal, but allows 
the spectroscopic interpretation based on Koopmans' theorem (the Hartree-Fock scheme). 
Another representation is local and can be directly applied to the analysis of interatomic 
magnetic interactions (the Kohn-Sham scheme). 



VII. SUMMARY 



Being based on the one-band extended Hubbard model, we have discussed several aspects 
of the Kohn-Sham formulation of the density-functional theory in solids. Although the con- 
clusions are strictly applicable only to the mean-filed solutions of this particular model, the 
analysis presented in the paper suggests that the old-standing dilemma between the descrip- 
tion of the ground-state properties in the density-functional theory and the description of (at 
least some) one-particle excitations may be resolved in a more rational way. Particularly, we 
have shown that the conductivity gap and the interatomic magnetic interactions relevant to 
the spin-wave excitations for this model are entirely defined by the parameters of the ground 
state. This therefore means that in this case there is no fundamental difference between the 
ground-state properties, the conductivity gap and the spin-wave excitation energies, and 
all of them can be described by the same set of one-electron Kohn-Sham equations. The 
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KS eigenvalues may not be necessarily coincide with the true one-particle excitation ener- 
gies, depending on the form in which the nonlocal exchange interactions contribute to each 
particular type of excitations. 

Thus, we hope that results of our model analysis may provide a new insight into the 
interpretation of one-electron equations underlying the Kohn-Sham scheme and stimulate 
new researches in this direction. As an extension of the present analysis, we believe that 
it would be particularly interesting: (i) to go beyond the mean-field approximation and to 
take into account the correlation effects; (ii) to go beyond the one-band picture, to take into 
account the effects of orbital degeneracy and to clarify the roles played by intersite Coulomb 
interactions in the problem of orbital polarization in solids.— 
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TABLE I: Symmetry constraints superimposed on the density matrix and the KS potential. The 
uniform potential shift is set as zero. Note also that the number of electrons is ne=noTv[p^). In 
the COF state, the minority-spin bands are requested to be unoccupied; in all other respects the 
parameters S j 2 can be arbitrary. 

state density matrix KS potential 



AF p 
CON 
COF p^=0 



d^pa^ -BI=b{=B. 
T 



-°2 



-bI=bI 



f{k)=cosk 



f{k)=cosk+cosk 




f(k)=cosk +cosk +cosk 

J \ ' X V z 



f{k)=4cos{k/2)cos{k /2)cos{k/2) 




FIG. 1: Doubling of the unit cell associated with the spin or/and charge ordering. Two nonequiv- 
alent sites {u=l,2) are shown by white and black circles. Within a doubled unit cell, the site 1 is 
located in the origin and the position of the site 2 is defined by the vector r. z is the coordination 
number, a is the lattice parameter, /(k) is the structure factor - see text (k is in the units of 1/a). 
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FIG. 2: Tx versus for various lattices and ne=no. 
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FIG. 3: Characteristic behavior of Tx, Tz and b calculated self-consistently as the function of the 



electron number near ne=no. 
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